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Within a hydrodynamical approach wc investigate the sensitivity of single inclusive mo- 
mentum spectra of hadrons in 158 A GeV Pb+Pb collisions to three different equations of 
state of nuclear matter. Two of the equations of state are based on lattice QCD results 
and include a phase transition to a quark-gluon plasma. The third equation of state has 
been extracted from the microscopic transport code RQMD under the assumption of com- 
plete local thermalization. All three equations of state provide reasonable fits to data taken 
by the NA44 and NA49 Collaborations. The initial conditions before the evolution of the 
fireballs and the space-time evolution pictures differ dramatically for the three equations of 
state when the same freeze-out temperature is used in all calculations. However, the softest 
of the equations of state results in transverse mass spectra that are too steep in the central 
rapidity region. We conclude that the transverse particle momenta are determined by the 
effective softness of the equation of state during the fireball expansion. 

PACS numbers: 24.10.Jv, 24.85.+p, 25.75.+r 



> 

o 
o 

OS 



o 

=5 



X 



A goal of studying ultrarelativistic heavy ion collisions is 
to understand nuclear matter under extreme conditions 
such as high energy density, high particle density, and 
whether a transition from a quark-confined hadron state 
to the deconfined state (quark-gluon plasma, QGP) takes 
place (see, e.g., ref. jl]]). A highly excited state of mat- 
ter and a possible transition to a QGP will occur at the 
early stages of the heavy-ion collisions; these differences 
will be reflected in different equations of state for nuclear 
matter. These in turn will affect the maximum energy 
density and compression achieved in a heavy ion reac- 
tion. To investigate the nature of the hadronic equation 
of state and possible phase transitions, it is very impor- 
tant to obtain information about the conditions that ex- 
ist at maximum compression, the subsequent evolution 
of the system, and the final free hadrons that can be 
measured in the detectors. 

The equation of state, (EOS), is part of a hydrodynamic 
model, and the connection between the EOS and exper- 
imental observables is relatively direct. In order to ad- 
dress issues rased above, we shall use the simulation code 
HYLANDER-C (c/., ref. J| and refs. therein), with dif- 
ferent equations of state and Cooper- Frye freeze-out ||, 
to predict the single particle momentum distributions. 
The experimental spectra from the 158A GeV Pb+Pb 
collisions, measured by the NA44 j|j5) and NA49 Collab- 
orations H , are used as fit criteria for the hydrodynam- 
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ical (one-fluid type) calculations. 

In our approach, virtually any type of EOS can be consid- 
ered by solving the relativistic Euler-equations |Q. The 
coupled system of partial differential equations necessary 
that describe the dynamics of a relativistic fluid (with 
given initial conditions) contains an equation of state 
which we write in the form 



P(e, n) = c 2 (e, n) e . 



(1) 



In eq. ([j]), the quantities P, e, and n are the pressure, 
the energy density, and the baryon density, respectively. 
The proportionality constant, c 2 , is in general a function 
of e and n. In the following, we shall assume that in 
particular the n dependence is negligible for the energy 
regime under consideration. From only the knowledge of 
the speed of sound, c 2 (e), one can then calculate c 2 (e) by 
evaluating the integral S 



c 2 (e) 



1 



cl(e>)de> 
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Assuming an adiabatic expansion, the temperature, T(e), 
can be calculated from I01 



T(e) = T exp 



cl(e>) de> 



„ (1 + c 2 (e')K 



(3) 



where T = T(eg) has to be specified for an arbitrary 
value of cq. 

Three equations of state are discussed in this study. The 
first one, EOS-I, which we use in the following exhibits 
a phase transition to a quark-gluon plasma at a critical 
temperature T c = 200 MeV with a critical energy density 
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e c = 3.2 GeV/fm 3 (c/., ref. §, and the rcfs. in §). The 
second equation of state, EOS-II, is also a lattice QCD- 
based EOS which has recently become very popular in the 
field of relativistic heavy-ion physics (c/., ref. juj). This 
equation of state includes a phase transition at T c = 160 
MeV with a critical energy density e c w 1.5 GeV/fm 3 . 
The third equation of state, EOS-III, has been extracted 
from the microscopic transport model RQMD |ll]] un- 
der the assumption of complete local thermalization [jl2] , 
and does not include a transition to a QGP. In Fig. 1 
the three equations of state are plotted in many differ- 
ent representations. For instance, in Fig. 1(d) the plot 
of c 2 = P(e)/e emphasizes the existence of a minimum 
in P/e at e = e c = 3.2 GeV/fm 3 (w 1.5 GeV/fm 3 ) for 
EOS-I (EOS-II), referred to as the softest point of the 
EOS. It corresponds to the boundary between the gener- 
alized mixed phase and the QGP JhJ. As seen in Fig. 1, 
EOS-II yields a much softer equation of state than does 
EOS-I. 

The parametrizations for the speed of sound which we 
have used here are |n| 



c 2 (e) = [a + ft tanh( 7 ln(e/e ) + 5)} 
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(i-£)r 2 



ln^(e/e ) + r2 



(4) 



in the case of EOS-I; 

c 2 D (e) = a + f3 tanh[7ln(e/eo)] 
-8 tanh[£hi(e/e')] , 

in the case of EOS-II; and 



c 2 (e) = i a- p (e/e ) 



e < e' 
e>e' 



(5) 
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in the case of EOS-III. The parameters of eqs. (^) - (||) 
and the choices for To = T(cq) are listed in Table I. 
In the following it is assumed that due to an experimental 
uncertainty for the centrality of the collision, only 90% of 
the total available energy and the total baryon number 
have been observed. It is then possible to find initial 
distributions for the three equations of state, such that 
one can reproduce the single inclusive momentum spectra 
of 158 AGeV Pb+Pb collisions (c/., Fig. 3). 
Fig. 2 shows the initial distributions for the energy den- 
sity, e(z), the baryon density, ub(z), and the fluid rapid- 
ity, ijf(z), plotted against the longitudinal coordinate z. 
We use here for the initial distributions the initial condi- 
tion scenarios which have been extensively described in 
refs. |l3| and [0|. In particular, it is assumed that an 
initial transverse fluid velocity field is absent, and the ini- 
tial longitudinal distributions for energy density, e, and 
baryon density, hb, are smeared out with a Woods-Saxon 
parametrization in the transverse direction, rj_ (c/., ref. 
fl5l|). The choices for the initial parameters, which in 



each case provide the best fit results for the hadronic 
single inclusive momentum spectra of the two considered 
experiments, are listed in Table II. 

The maximum initial energy density is 6a = 15.3 (6.55 
(8.66)) GeV/fm 3 , and the maximum initial baryon den- 
sity is n B lax = 3.93 (1.24 (1.92)) fm~ 3 , for EOS-I (EOS- 
II (EOS-III)), respectively; 71% (30% (36%)) of the bary- 
onic matter is initially located in the central fireball re- 
gion. Hence, EOS-I predicts a much larger stopping than 
EOS-III does, and EOS-III predicts a much larger stop- 
ping than EOS-II does. 

Fig. 3 shows results of the hydrodynamical calculations 
compared to the single inclusive momentum spectra of 
the 158 AGeV Pb+Pb collisions which have been mea- 
sured by the NA44 {§, f§ and NA49 Collaborations [§. 
All single inclusive momentum spectra have been evalu- 
ated in the nucleus-nucleus center of mass system (y cm 
= 2.91), and for all calculations the freeze-out tempera- 
ture Tf = 139 MeV has been used. The energy density 
at freeze-out is e f = 0.292 (0.126 (0.130)) GeV/fm 3 for 
EOS-I (EOS-II (EOS-III)). In particular, the freeze-out 
energy density, e/, is a factor ~ 2.2 larger in the calcu- 
lation using EOS-I compared to the other two calcula- 
tions. Inspecting Fig. 4, we found that the total lifetime 
of the fireball is reduced by approximately this factor 
in the calculation using EOS-I. Fig. 4 also shows that 
EOS-II and EOS-III yield fireballs of similar lifetimes 
(cf. also Table II). These particular features show up 
in Bose-Einstein correlation (BEC) calculations (cf. ref. 
p|); more on BEC will be published elsewhere Jig . The 
relatively strong effect in two-particle correlation func- 
tions was also observed in a recent report |l7j . 
All three calculations yield resonable agreement with 
both experiments. However, the calculations for the very 
soft EOS-II yield transverse mass spectra that are too 
steep in the central rapidity region. The magnitudes of 
the slopes in the transverse mass spectra have their origin 
in the transverse velocity fields at freeze-out. For EOS-II 
its corresponding transverse velocity field is the weakest 
(cf. Table II) of the three EOS. In fact, it is about 30% 
smaller than the transverse velocity fields for EOS-I and 
EOS-III. The final velocity is proportional to the time 
integral 



/ P 



(7) 



-(f) dt' , 



where Pj e is proportional to the fluid acceleration |18| , |l9| . 
In Fig. 1(d) we exhibit for each EOS we used the starting 
values of P / e with respect to the initial maximum energy 
density £a at transverse position r± = 0, and the final 
values of P/e at breakup energy densities, e/. For the 
given total fireball lifetime, t max = tf — tt, the effective 
P/e of EOS-II is obviously too small, whereas P/e of 
EOS-I and EOS-III appear to have the correct effective 
softness. We stress that the total lifetime, t max , in the 
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case of EOS-I is reduced by a factor ~ 2.1 compared to 
t max of EOS-III, but the ratio P/e of EOS-I is effectively 
increased for the considered expansion period. 
In summary, we have tested three types of equation of 
state for 158A GeV Pb+Pb collisions. We found that 
the measured single inclusive momentum spectra are sen- 
sitive to both the initial conditions and to the equation 
of state. EOS-I and EOS-III both provide reasonable fits 
to the data. For the EOS-II a softest point is assumed 
and its corresponding velocity field is the weakest. As a 
result, the transverse mass distributions calculated from 
this equation of state are found to be too steep in the cen- 
tral rapidity region as compared to the experimental ob- 
servations. Obviously, the transverse particle momenta 
are sensitive to the effective softness of the equation of 
state during the fireball expansion. 

As mentioned earlier, our conclusions are based on one- 
fluid hydrodynamic model calculations. Furthermore, 
the freeze-out temperature is fixed at 139 MeV for all 
three cases under study. As long as we limit our compar- 
ison of different scenarios within the same framework, our 
conclusions are valid. Still, several crucial questions like 
whe re and when the collective velocity field [f20| is de- 
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velo ped during the collision, and whether it develops in 
the QGP phase or in the hadronic phase are not answered 
yet. In addition, the exact condition for the system evolv- 
ing from interacting hadrons to freeze-streaming ones is 
not known. To shed light in this direction, an even more 
sophisticated analysis is needed. We believe that the 
combined analysis of single inclusive momentum spectra 
and Bose-Einstein correlations will help to further con- 
strain |2lj the limits of possible candidates for the equa- 
tion of state of nuclear matter in relativistic heavy-ion 
collisions. 
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FIG. 1. Energy density, e, ratio of pressure and energy den- 
sity, P/e, speed of sound, Cq, and temperature, T, as func- 
tions of T and/or e, for the equations of state EOS-I (solid 
lines), EOS-II (dashed lines), and EOS-III (dotted lines), re- 
spectively (see text). The open circles in plot (d) correspond 
for each EOS to the starting values of P/e with respect to 
the achieved initial maximum energy density ea at transverse 
position r± = 0, whereas the filled squares represent the final 
values of P/e at breakup energy densities, ej. 
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FIG. 2. Initial distributions of (a) the energy density, e, 
(b) the baryon density, ns, and (c) the fluid rapidity, yp, 
plotted against the longitudinal coordinate z at transverse 
position r± — 0. The solid (dashed (dotted)) lines indicate the 
initial distributions for the calculation using EOS-I (EOS-II 
(EOS-III)). 
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FIG. 3. (a) Rapidity spectra and (b) transverse mass spec- 
tra, l/m±dN/dm±dy, of negative hadrons, bT , (c) rapidity 
spectra of protons (without contributions from A decay) and 
(d) transverse mass spectra, l/m±dN/dm±, of protons (in- 
cluding contributions from A decay), p, positive kaons, K + , 
and positive pions, ir + , respectively. The solid (dashed (dot- 
ted)) lines indicate the results of the calculations when using 
equation of state EOS-I (EOS-II (EOS-III)). The open circles 
represent preliminary data taken by the NA49 Collaboration, 
whereas the filled squares represent final data taken by the 
NA44 Collaboration. 
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FIG. 4. Isothermes for the relativistic fluids governed by 
EOS-I, EOS-II and EOS-III at r± = 0, respectively. The 
outer lines correspond to a temperature, T = 139 MeV, and 
each successively smaller curve represents a reduction in tem- 
perature by AT = 20 MeV. The lines c = 1 represent the 
light cone. 

TABLE I. EOS Parameters. 



EOS-I 


EOS-II 




EOS-III 


a 4/15 


0.2582 




0.15 


P 1/15 


1/6 




0.02 


7 0.24 


10.0 




— 


8 1.05 


0.0915 




2.00 


£ 0.15 


1.50 




- 


T 2 0.73 






- 


e [GeV/fm 3 ] 


0.27 




2.80 


e [GeV/fm 3 ] 3.20 


1.45 




1.00 


T [MeV] 200 


160 




180 


TABLE II. Properties 


of the fireballs. 






EOS-I EOS-II 


EOS-III 


Initial parameters 






Rel. fraction of thermal energy 


0.55 


0.20 


0.28 


in the central fireball, Kl 








Longitudinal extension of the 


1.20 


1.00 


1.10 


central fireball, A [fm] 








Rapidity at the edge of the cen- 


1.00 


0.85 


0.90 


tral fireball, j/a 








Rapidty at maximum of initial 


0.80 


1.50 


1.10 


baryon y distribution, y m 








Width of initial baryon y distri- 


0.32 


1.00 


0.60 


bution, a 








Output 








Max. initial energy density, ea 


15.3 


6.55 


8.66 


[GeV/fm 3 ] 








Max. initial baryon density, 


3.93 


1.24 


1.92 


nT x [fm~ 3 ] 








Rel. fraction of baryons in cen- 


0.71 


0.30 


0.36 


tral fireball, f„ B 








Freeze-out energy density, e/ 


0.292 


0.126 


0.130 


[GeV/fm 3 ] 








Max. transverse velocity at 


0.46 


0.30 


0.42 


freeze-out, v^ ax [c] 








Lifetime of fireball, t ma x [fm/c] 


13.1 


29.3 


27.3 


Lifetime of QGP, tQap [fm/c] 


2.4 


3.0 
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